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The Born-Oppenheimer approximation (BO) [lj is the standard 
ansatz to describe the interaction between electrons and nuclei. 
BO assumes that the lighter electrons adjust adiabatically to the 
motion of the heavier nuclei, remaining at any time in their in- 
stantaneous ground-state. BO is well justified when the energy 
gap between ground and excited electronic states is larger than the 
energy scale of the nuclear motion. In metals, the gap is zero and 
phenomena beyond BO (such as phonon-mediated superconductiv- 
ity or phonon-induced renormalization of the electronic properties) 
occur [2]. The use of BO to describe lattice motion in metals is, 
therefore, questionable pSH]. In spite of this, BO has proven effec- 
tive for the accurate determination of chemical reactions |5j, molec- 
ular dynamics [6], [7] and phonon frequencies [9], [HI 110] in a wide 
range of metallic systems. Graphene, recently discovered in the 
free state |11L I12j. is a zero band-gap semiconductor [13J , which be- 
comes a metal if the Fermi energy is tuned applying a gate-voltage 
V g | |14L 112] . Graphene electrons near the Fermi energy have two- 
dimensional massless dispersions, described by Dirac cones. Here 
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we show that a change in V g induces a stiffening of the Raman 
G peak (i.e. the zone-center E 2s optical phonon |15L 116]). which 
cannot be described within BO. Indeed, the E 29 vibrations cause 
rigid oscillations of the Dirac-cones in the reciprocal space [17J. If 
the electrons followed adiabatically the Dirac-cone oscillations, no 
change in the phonon frequency would be observed. Instead, since 
the electron-momentum relaxation near the Fermi level |18|, 1191 120] 
is much slower than the phonon motion, the electrons do not fol- 
low the Dirac-cone displacements. This invalidates BO and results 
in the observed phonon stiffening. This spectacular failure of BO 
is quite significant since BO has been the fundamental paradigm 
to determine crystal vibrations from the early days of quantum 
mechanics fTl [9l l2T| \8\ [TO] . 

Graphene samples are prepared by micromechanical cleavage of bulk 
graphite at the surface of an oxidized Si wafer with a 300 nm thick ox- 
ide layer, following the procedures described in Ref. [11]. This allows us 
to obtain graphene monocrystals exceeding 30 microns in size. By using 
photolithography, we then make Au/Cr electrical contacts, which enable 
the application of a gate voltage, V 9 , between the Si wafer and graphene 
(Fig. [T]A.,B). The resulting devices are characterized by electric-field-effect 
measurements [T2l [HJ [22], yielding a charge carrier mobility fi of 5,000 to 
10,000 cm 2 /Vs at 295K and a zero-bias (V 5 =0) doping of ~10 12 cm" 2 [23]. 
This is reflected in the existence of a finite gate voltage V n at which the 
Hall resistance is zero and the longitudinal resistivity reaches its maximum. 
Accordingly, a positive (negative) V 9 -V n induces an electron (hole) doping, 
having an excess-electron surface-concentration of n=r/(V g - V n ). The coeffi- 
cient rj &7.2 10 10 cm~ 2 /V is found from Hall effect measurements and agrees 
with the geometry of the resulting capacitor [T2l dU fH"] . 

Unpolarized Raman spectra are measured at 295 and 200 K in ambient air 
and in vacuum (<5 10 -6 mbar), respectively, with a Renishaw spectrometer 
at 514nm using a 50 x long working distance objective, Fig. [TJ3. The incident 
power is kept well below 4mW in order to avoid sample damage or laser 
induced heating [T5]. The Raman spectra are measured as a function of the 
applied V 9 , Fig. [2A. Each spectrum is collected for 30 seconds. The applied 
gate voltage tends to move the Dirac point, especially at room temperature. 
We thus determine the V g corresponding to the minimum G peak position, 
and use this to estimate V n . Fig. [2]A.,B show that the G peak upshifts with 
positive applied V 9 -V n . It also shows a similar trend, albeit over a smaller 
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Figure 1: (A) Optical micrograph of the contacted graphene sample. (B) 
Schematic of the Raman+transport measurement system 

voltage range, for negative V 9 -V n . This upshift for both electron and hole 
doping is qualitatively similar to that reported by Yan et al. for electrically 
doped graphene measured at 10K [24]. 

The Raman G peak of graphene corresponds to the E 2g phonon at T [T5} 
IT6] . Phonon calculations done within BO for undoped graphene and graphite, 
show the presence of a Kohn anomaly in the phonon dispersion of the E2 9 
mode near T [25]. A Kohn anomaly is the softening of a phonon of wavevector 
q ~ 2kp, where kp is a Fermi surface wavevector [26]. By doping graphene, 
intuitively one could expect that the change in the Fermi surface should 
move the Kohn anomaly away from q=0 and, thus, stiffen the T phonon 
detected by Raman measurements, which would be in agreement with our 
experiments. To validate this picture, we need to compute the frequency of 
the E 2g mode in doped graphene. 

In graphene, the electronic bands near the high-symmetry K points are 
well described by a Dirac dispersion [T3] e(k, n*) = Tivpk and e(k, n) = 
—hvpk, where k + K is the momentum of the Dirac Fermions, vf is the 
Fermi velocity and Hvf = 5.52 eV A, from density functional theory (DFT) 
calculations [25] (Fig. [3]A). The Dirac point is defined by the crossing of 



3 



1580 1585 1590 1595 -6 -4 -2 2 4 6 
Raman Shift (cm 1 ) Electron Concentration (10 12 /cm 2 ) 

16 
15 
14 
13 
g 12 

I 11 

Li. 

10 

9 
8 
7 

-6 -5 -4 -3 -2 -1 1 2 3 4 5 6 ~" -6 -5 -4 -3 -2 -1 1 2 3 4 5 6 
Electron Concentration (10 12 /cm 2 ) Electron Concentration (10 12 /cm 2 ) 

Figure 2: (A) G peak of graphene measured at room temperature as a func- 
tion of V g . The red spectrum corresponds to the Dirac point. (B-C) G peak 
position as a function of electron concentration at 200 and 295 K: (black 
dots) measurements; (red-dashed line) adiabatic Born-Oppenheimer; (blue 
line) finite-temperature non-adiabatic calculation from Eq. [6j (thin dashed 
black) simplified non-adiabatic calculation from Eq. [5j (D) G peak Full 
Width at Half Maximum, FWHM(G), at 200K as a function of electron con- 
centration: (dots) measured; (Blue line) theoretical FWHM of a Voigt profile 
obtained from a Lorentzian component given by Eq. and a constant Gaus- 
sian component of ~8 cm -1 . 
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these conic bands and coincides with K, Fig. [3]A.. Thus, at zero temperature, 
the doping-induced shift of the Fermi level from the Dirac-point is e F = 
sgn(n)y/mrhv F , where sgn(x) is the sign of x. 

The E 2g phonon in graphene consists of an in-plane displacement of the 
carbon atoms by a vector ±u/\/2 as sketched in Fig. 03. In presence of such 
atomic displacements, the bands are still described by a cone (i.e. a gap does 
not open) with the Dirac-point shifted from K by a vector s (Fig. GB,C) [17J. 
In practice, the atomic-pattern of the E 2g vibrations is mirrored into an 
identical pattern of Dirac-point vibrations in the reciprocal space. The de- 
pendence of the electronic-bands on u can be obtained (see supplementary 
information) from the DFT electron-phonon coupling matrix-elements (Eq. 6 
and note 24 of Ref. [25]): 



e(k, ttVtt, u) = ±hv F \k-s(u)\ (1) 



where s ■ u = 0, s = u^2(Dl) F /(hv F ), and (L> 2 ) F = 45.6 (eV) 2 /!" 2 is the 
deformation potential of the E 2g mode [27]. Eq. [I] well reproduces the modi- 
fication of the DFT band structure of graphene due to a static displacement 
(frozen-phonon) of the atoms according to the G phonon pattern. 

The knowledge of the electronic-bands (in the presence of a phonon) al- 
lows the determination of the phonon energy hu €F as a function of e F . In 
particular, 

. . . t h d 2 AE 

= ^ " ^ = 2Mu (duf ' (2) 

where M is the carbon mass, Auj <C uj and AE is the variation of the 
electronic energy with e F . 

Within BO, AE{u) is computed assuming a static atomic displacement. 
Under this hypothesis, for any given displacement u, the electrons are sup- 
posed to be in the ground state configuration, i.e. the bands e(k, 7r*,u) are 
filled up to e F (Fig. EB). Thus, the adiabatic AE is 

44 r 

AE(u) = -— e(k,vr*,u)d 2 fc, (3) 

{ZTl z ) Je(k,vr*,u)<e F 

where we consider e F > 0, A = 5.24 A 2 is the unit-cell factor 4 

accounts for spin and K-point degeneracy. Combining Eq. [1] and El we 
have that AE does not depend on u and TiAu = 0. Thus, within BO, 
the Raman G peak position is independent of e F , in contrast with the trend 
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Figure 3: Schematic it band structure of doped graphene near the high sym- 
metry K point of the Brillouin zone. The filled electronic states are colored 
in green. (A) Bands of the perfect crystal. The Dirac point coincides with K, 
the electronic states are filled up to the Fermi energy ep and the Fermi surface 
is a circle centered at K. (B) Bands in presence of an E2 5 lattice distortion. 
The Dirac points are displaced from K by ±s. In the adiabatic approxima- 
tion, the electrons remain in the instantaneous ground state: the bands are 
filled up to ep and the Fermi surface follows the Dirac-point displacement. 
The total electron-energy does not depend on s. (C) Bands in presence of 
an E 2g lattice distortion. In the non-adiabatic case, the electrons do not 
have time to relax their momentum (through impurity, electron-electron and 
electron-phonon scattering) to follow the instantaneous ground state. In ab- 
sence of scattering, the electron momentum is conserved and a state with 
momentum k is occupied if the state with the same k is occupied in the 
unperturbed case. As a consequence, the Fermi surface is the same as in 
the unperturbed case and does not follow the Dirac-cone displacement. The 
total electron-energy increases with s 2 resulting in the observed E2 9 -phonon 
stiffening. (D) Atomic pattern of the E 29 lattice distortion, corresponding to 
the Raman G peak. The atoms are displaced from the equilibrium positions 
by ±u/^/2. Note that the displacement pattern of the Dirac points (in re- 
ciprocal space) is identical to the displacement pattern of the carbon atoms 
(in real space). 6 



reported in Fig. 06, C. Note that Ref. |28J calculated the doping dependence 
of the G peak position within BO, but found significant phonon softening 
for increasing doping, quite the opposite of the experiments, and in contrast 
with our BO calculation. 

The failure of the frozen-phonon calculation, urges us to re-examine the 
assumptions underlying BO. The E 2g phonon is a dynamical perturbation 
described by a time-dependent lattice displacement u(t) = ucos(a>ot) oscil- 
lating at the G peak frequency. Within BO, it is assumed that, at any given 
time t, the electrons are in the adiabatic ground state of the instantaneous 
band structure e(k, n*, u(t)). However, the period of the G peak vibrations is 
~21 fs, which is much shorter than the typical electron- momentum relaxation 
times r m (due to impurity, electron-electron and electron-phonon scattering). 
Indeed, a r m of a few hundreds fs is deduced from the electron-mobility in 
graphene [18] and from ultra-fast spectroscopy in graphite [19], [20] . As a 
consequence the electrons do not have time to relax their momenta to reach 
the instantaneous adiabatic ground state, as assumed in BO. The departure 
from the adiabatic ground state can be accounted for in the calculation of 
AE, by filling the perturbed bands, e(k, 7r*,u) with the occupations of the 
unperturbed bands e(k, ir*,0), as in Fig. [3D: 

AA r 

AE(u) = -—- / e(k, vr*, u) d 2 k + 0(u 3 ). (4) 

This equation is valid in the limit e F ^> huo/2, and can be rigorously derived 
using time dependent perturbation theory (see supplementary information). 
In this case, the non-adiabatic energy, AE, depends on u. By replacing Eq. [TJ 
and Eq. H]in Eq. [2] and performing the integral we get: 

TiA(D 2 T ) F . . „ . 

hALJ = M ft V 6F = a 6F ' 5 

nMuj (hv F ) z 

where a' = 4.39 KT 3 . 

The result of Eq. [5] can be extended to any e F and finite temperature T 
using time dependent perturbation theory [29] to obtain: 

p, a r [/(e - ^f) - /(e)]e 2 sgn(e) 

hAuj = a P / — — -y— de, (6) 

where P is the principal part, and / is the Fermi-Dirac distribution at T. 
Fig. [2j3,C shows the excellent agreement of our non-adiabatic finite T calcu- 
lation (Eq. [6]) with the experiments. The measured trends are also captured 
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by the simplified model, Eq. By comparing the predictions of the BO 
calculation and of the non-adiabatic model, we conclude that the stiffening 
of the E 2g mode with \ep\ is due to the departure of the electron population 
from the adiabatic ground state. 

A pictorial interpretation of this phenomenon can be obtained by con- 
sidering what happens to a filled glass when shaken horizontally The liquid 
gravitational-energy and its level mimic the electronic energy AE and ep, 
respectively. The shaking frequency mimics the phonon frequency and the 
relaxation time of the liquid-surface mimics the electron relaxation time. 
If the motion of the glass is slow, the liquid surface remains flat and its 
gravitational-energy is independent of the glass horizontal position, as in 
Eq. |3] and in Fig. [3B. If the motion of the glass is rapid, the liquid surface 
profile is not flat and its gravitational-energy increases with the displacement 
of the glass, as in Eq. 0] and Fig. [3D. To push the analogy even further, one 
should use a non-cylindrical glass, for which the liquid surface increases with 
the liquid level. In this case, the higher the liquid-level, the larger the dif- 
ference between the gravitational energies in the fast- and slow-shaken glass. 
Indeed, in graphene, the higher the Fermi level, the larger the difference be- 
tween the non-adiabatic AE and the adiabatic AE. This causes the observed 
stiffening of the phonon frequency with ep. 

The validity of our model is further confirmed by the analysis of the 
G peak linewidth. The phonon decaying into an electron-hole pair gives 
the most important contribution to the homogeneous broadening of the E 2g 
phonon. The full-width at half-maximum, 7, can be computed extending to 
finite T and ep 7^ the results of Ref. [27] : 



where c is the speed of light. At T = 0, 7 = 11 cm -1 for = and 7 
drops to zero for ep > Huiq/2 because the scattering process is forbidden by 
the Pauli exclusion principle [27] . Fig. [2p shows a good agreement between 
the experimental and theoretical 7, once a constant inhomogeneous Gaussian 
broadening of ~ 8 cm -1 is added to the electron-phonon contribution of Eq. [7J 
Concluding, the observed stiffening of the E 2g phonon in doped graphene 
represent a spectacular failure of the adiabatic Born-Oppenheimer approx- 
imation. Within BO, the energy of a zone-center phonon is determined by 
two contributions: the distortion of the electronic bands, associated with the 
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phonon displacement, and the consequent rearrangement of the Fermi sur- 
face. These two contributions cancel out exactly in graphene because of the 
peculiar rigid motion of the Dirac-cones, associated to the E2 S phonon. In 
general, a correct phonon treatment should not include the BO Fermi-surface 
rearrangement, whenever the electron-momentum relaxation time (near ep) 
is longer than the phonon period. We anticipate the failure of BO, shown 
here, to have important consequences in the description of vibrational prop- 
erties of carbon-nanotubes and in phonon-mediated superconductors. 
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Supplementary Information 



Derivation of Eq. [T] 

The electronic Hamiltonian for the n,ir* basis can be written as a 2 x 2 matrix: 



H(k,0) 



fivpk 
—hvpk 



(8) 



where k is a small in plane wave- vector and K + k is the electronic momen- 
tum. Let us consider a distortion of the lattice according to a T — E 2g phonon 
pattern (note that the T — E 2g phonon is doubly degenerate). At the lowest 
order the 7r-bands Hamiltonian changes as 



H(k,u) = H(k,0) + dH ^ 0) u 

ou 



(9) 



where u is the phonon normal coordinate (the two atoms in the unit-cell are 
displaced by ±u/\^2 along a given direction in the plane). dH/(du) can be 
obtained from the ab-initio deformation potential matrix elements. Following 
Ref. [25] (Eq.6 and note 24) and Ref. [27], for the E 2g phonon mode and for 
a small k 



fkvr* 



OH 

du 



kvr*) 



(Dl) F [l + cos(26)] (10) 



M -ik.) 



(££) F [l-cos(20)], 



(11) 



where |k7r / 7r*) are the electronic states with momentum K+k and 9 is the 
angle between k and the direction perpendicular to the atomic vibration. 
Taking the square root of Eqs. [TO], [11] and inserting them into Eq. [9] 



H(k,u) 



Tivpk 






-hvpk 



2{Dl) 1 



cos(^) Bin(0) 
sin(^) — cos(^) 



e = ±hv F yk 2 + s 2 + 2ks cos(6 ) ) = ±fiv F \k — s(u) 
where s is defined in the main text. 



u. 



The eigenvalues of Eq. [12] are then 



;i2) 



(13) 
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Derivation of Eq. 31 

Considering the Taylor expansion of AE(u) in u, Eq. H]is equivalent to: 

d 2 A , s d 2 { AA f „ x ] , , 

A£( M ) = — — — — / e(k,7r*,u)d 2 £; . (14) 



(du) 2 1 (rfu) 2 [_ (27r) 2 7 e (k,7r*,0)<e F 

In this section, we will demonstrate that 

h d 2 ( 4A r x , 9 , 1 . . 

frAu = - - — — ^ - — — / e(k,7r*,u) d 2 k }, (15) 

at T = and under the condition 

t F > Uujq/2. (16) 

Using Eq. El Eq. [4] will then immediately follow. 

Within time dependent perturbation theory, TiAu is (see Eq. 10 of 
Ref. [29]): 

hAu = _5_[JS'(wb) - ifj(wb)]. (17) 
2Mu;o 

were at T = 0, 

= ^- £ lAco, k e| 2 { + ! — T- ) • (18) 

Here the index o and e denotes the occupied (eko < cf) and empty bands 
(eke > £f), and 

^ko,ke= (ko| — |ke). (19) 
aw 

Now we consider only the 7r and 7r* bands and we substitute 1/iVkXik with 
A/(27t) 2 J d 2 k, where A is the unit-cell area and the integral is restricted 
on a circle of radius k, centered on K. Assuming a Dirac dispersion for 
the 7r and 7r* bands, ek e — > 2e F . Thus, if the condition of Eq. [TBI holds. 
I^ke - e ko| huo and the Tvujq in the denominators of Eq. [18] can be neglected. 
Eq. [18] becomes 

F €F = 8A [ |-Pk7r»,k7r| 2 / 2Q\ 

(2vr) 2 i e (k,^,o)>e F ,fc<fc e(k,vr,0) -e(k,7r*,0) ' V ' 
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where e(k, ir/n*, 0) are the bands of the undistorted graphene structure. 
From Eq. [T7] 



n 



8A 



k7T 



2MlUq { (27r) 2 Je(k,n*,0)<e F e(k, 7T, 0) - e(k, 7T*, 0) 

From textbook static second order perturbation theory 

I n I 2 



eTfc 



1 d 2 e(k,7r*,u) 
2 



(d«) 2 

Substituting Eq. [22] in Eq. [21] we have 



u=Q e(k,7r*,0)-e(k,7r,0) 



n 



4A 



d 2 e{k,7c*,uj , 



2Mu { (2tt) 2 Je(k,n*,o)<e F (du) 2 



d 2 k 



(21) 



(22) 



(23) 



Eq. [15] is, finally, obtained by taking the derivation with respect to u in Eq. 
outside the integral. 
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